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Incremental  Diffraction  Coefficients  for  Planar  Surfaces,  Part  II: 
Calculation  of  the  Nonuniform  Current  Correction  to  PO 

Reflector  Antenna  Patterns 


1.  INTRODUCTION 

Numerous  computer  codes  exist  for  calculating  the  far  fields  of  reflector  antennas  by  integrating 
either  the  physical  optics  (PO)  current  or  the  aperture  fields  obtained  from  geometrical  optics  (GO).*'^ 
Both  the  PO  current  integration  and  the  GO  aperture-field  integration  produce  far  fields  with  reasonable 
and  comparable  accuracyt  out  to  a  sidelobe  angle  that  depends  upon  whether  the  edges  of  the  reflector 
are  sharp,  flared,  absorber-lined,  serrated,  etc.^  In  particular,  for  antennas  with  locally  wedge-shaped 
conducting  edges,  both  PO  current  and  CO  aperture-field  integration  predict  nearly  the  correct  copolar 
fields  out  to  20®  or  more.^  However,  even  for  perfectly  conducting  reflectors  with  wedge-shaped  edges, 
the  PO  current  or  GO  aperture-field  integrations  neglect  the  nonuniform  currents’  or  fields^  near  the 
edges  and  thus  fail  to  produce  reliable  cross-polarized  fields,  far-out  sidelobes  of  the  copolarized  fields 
and  fields  near  nullsJ  Thus  a  number  of  techniques  have  been  used  to  improve  the  accuracy  of  PO 
current  and  GO  aperture-field  integrations. 

The  geometrical  theory  of  diffraction*  (GTD)  or  the  physical  theory  of  diffraction’  (PTD)  can  be 
used  to  calculate  more  accurately  the  far-out  sidelobes  of  reflectors.  However,  GTD  introduces  the 
problem  of  matching  smoothly  the  far-out  sidelobes  to  the  near-in  sidelobes,  and  both  PTD  and  GTD, 
when  applied  directly  to  obtain  the  far-fields,  fail  (for  a  general  illumination)  to  produce  accurate  co- 
and  cross-polarized  fields  in  the  mainbeam  region  (which  is  a  focal  region  for  the  GTD  and  PTD 
diffraction  coefficients)  Moreover,  if  conventional  GTD  or  PTD  is  used  to  first  calculate  fields  on  a 
near-field  surfree  before  performing  an  aperture-field  integration  to  obtain  the  far  fields,^-^  one  is  still 

(Received  for  publication  27  October  1987) 

fThe  accuracy  of  the  PO  current  and  GO  aperture-field  integrations  is  comparable,  provided  the  aperture  surfiice 
used  in  the  field  integration  caps  the  reflector.’ 
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foced  with  the  problem  of  computing  accurate  near  fields  in  the  focal  regions.  Also,  conventional  GTD  or 
PTD  is  not  well  suited  for  computing  end-point  contributions  such  as  from  corners,  even  if  higher  order 
distortions  of  the  current  near  the  end-points  are  neglected. 

To  overcome  this  problem  of  computing  helds  in  caustic  and  end-point  regions,  Mitzner*® 
introduced  PTD  “incremental  length  diffraction  coefficients,"  which,  when  multiplied  by  the  incident 
field,  could  be  integrated  along  the  rims  of  reflectors  (for  example)  to  obtain  the  diffracted  fields  of  the 
nonunifbrm  current  for  arbitrary  angles  of  incidence  and  scattering.  Mitzner  obtained  explicit  expressions 
for  the  incremental  diffraction  coefficients  of  the  perfeedy  conducting  wedge.  More  recendy,  Michaeli,*' 
unaware  of  Mitzner’s  report,  derived  expressions  for  the  GTD  incremental  diffraction  coefficients  of  the 
wedge,  and  Knott*^  has  shown  that  if  the  PTD  incremental  diffraction  coefficients  for  the  wedge  are 
subtracted  from  the  GTD  coefficients,  the  result  is  the  PO  incremental  diffraction  coefficients.  (The 
concept  of  incremental  diffraction  coefficients  was  used  implicitly  in  work  that  preceded  that  of 
Mitzner,’®  notably  that  of  Ufimtsev,*  Braunbek,^  Millar,’^"’^  {),£  “equivalent  (edge)  current”  work  of 
Ryan  and  Peters.’^  However,  this  previous  work  used  approximate  ad  hoc  expressions  for  the  incremen¬ 
tal  fields  that  did  not  necessarily  agree  with  the  diffracted  fields  outside  the  focal  regions.) 

Part  of  this  work  reviews  more  thoroughly  the  key  references  leading  to  GTD  and  PTD 
incremental  diffraction  coefficients,  and  provides  a  general  and  convenient  method  for  determining 
incremental  length  diffraction  coefficients  directly  from  the  conventional  two-dimensional  diffraction 
coefficients.  The  main  purpose  of  this  sequel  (Part  II)  is  to  provide  the  detailed  analysis  of  how  the  PTD 
incremental  diffraction  coefficients  for  the  half-plane  are  integrated  around  the  rims  of  reflectors  to 
obtain  the  far  fields  produced  by  the  nonuniform  currents.  These  computed  far  fields  of  the  nonuniform 
current  are  added  to  the  far  fields  computed  from  the  PO  current.  Excellent  agreement  with  the  far 
fields  obtained  from  a  method  of  moments  solution  to  the  electric  field  integral  equation  applied  to  a 
20-wavelength-diameter  reflector  with  an  8-wavelength  focal  length^-^®  shows  that  indeed  the  cross¬ 
polarization,  further-out  sidelobes,  and  fields  near  nulls  of  reflector  antennas  can  be  appreciably  changed 
by  the  fields  of  the  nonuniform  currents. 

Finally,  we  evaluate  the  fields  of  the  nonuniform  current  of  the  reflector  through  conventional 
asymptotic  evaluation  of  the  diffraction  integrals  and  compare  these  approximate  resuults  with  the  far 
fields  obtained  by  numerically  integrating  the  nonuniform  incremental  diffraction  coefficients. 


2.  ANALYSIS 

We  wish  to  improve  upon  the  PO  approximation  to  the  far  field  of  a  parabolic  reflector  antenna  by 
taking  into  account  the  fact  that,  in  the  close  vicinity  of  the  rim  of  the  reflector,  the  PO  current  does  not 
accurately  represent  the  total  current  excited  on  the  reflector  surface.  The  total  current  on  the  reflector 
surface,  K,  can  be  represented  as  the  sum  of  the  PO  current,  K'’°,  and  the  nonuniform  current,  K"“, 
which  is  prominent  only  near  the  reflector  rim  and  which  decays  rapidly  away  from  the  rim.  Correspond¬ 
ing  to  the  representation  of  the  total  current  as 


K  =  RPO  -h  K'"' 


is  the  representation  of  the  total  electric  far  field  scattered  by  the  reflector,  E*,  as  the  sum  of  the  PO  far 
field,  E**®,  and  the  far  field  of  the  nonuniform  current,  E"'*;  that  is. 
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(1) 


Ej  =  +  E^ 


with  the  far  field  given  in  terms  of  surface  integrals  of  the  respective  currents.  The  far-field  contribution  of 
the  nonuniform  current  can  be  found  by  subdividing  the  rim  into  differential  lengths,  determining  the  far 
field  of  each  such  differential  length,  and  integrating  this  far  field  over  the  rim  of  the  reflector. 

The  far  field  of  a  differential  length  (the  incremental  diffraction  coefficient)  can  often  be  well 
approximated  by  the  corresponding  far  field  of  a  canonical  scatterer  that  conforms  to  the  reflector  at  the 
rim  and  that  is  excited  by  a  plane  wave  whose  propagation  vector  is  in  the  direction  from  the  feed  to  the 
rim  element  and  whose  amplitude  and  phase  at  the  point  of  incidence  are  those  of  the  feed  radiation.  In 
general,  the  rim  of  the  reflector  may  not  conform  to  a  canonical  scatterer  with  a  known  incremental 
diffraction  coefficient.  However,  the  rims  of  many  reflectors  are  well  approximated  locally  by  perfectly 
conducting  wedges.  Incremental  diffraction  coefficients  for  the  perfectly  conducting  wedge  have  been 
derived  by  Mitzner'®  for  the  nonuniform  current,  Michaeli^^  for  the  total  current,  and  Knott^^  for  the  PO 
current.  (Since  the  total  current  is  the  sum  of  the  PO  and  nonuniform  currents,  any  two  of  the  coefficients 
determine  the  third.)  For  simplicity  here,  we  will  consider  the  special,  though  important,  case  of  a 
parabolic  reflector  whose  rim  can  be  modeled  locally  by  an  infinitesimally  thin,  perfectly  conducting 
half-plane.  The  correction,  to  the  PO  reflector  far  field  is  thus  obt^  ned  by  integrating  the  nonuniform 
current  incremental  diffraction  coefficient  for  the  half-plane  around  the  reflector  rim. 


2.1  Incremental  Diffraction  Coefficients  for  the  Nonuniform  Current  of  the  Half-Plane 
The  half-plane  nonuniform  current  incremental  diffraction  coefficient,  dE"“  (r),  is  given  by 


dEnu  (r)  =  dEs  (r)  -  dE^o  (fi 


(2) 


where  dE,  (r)  and  dE  (f)  are  the  half-plane  incremental  diffraction  coefficients  for  the  total  current  and 
PO  current  respectively,  and  r  is  the  vector  from  the  rim  point  to  the  far  field  point.  Expressions  for  dEs 
and  dE (r)  for  the  perfectly  conducting  wedge  were  derived  in  Part  I  and  shown  to  agree  with  Michaeli 
and  Knott  (and  thus  indirectly  with  Mitzner)  as  a  check  on  our  general  method  for  obtaining  incremental 
diffraction  coefficients.  (The  total  incremental  diffraction  coefficients  for  a  perfectly  conducting  half-plane 
were  also  derived  directly  in  Part  I.)  Specifically,  the  half-plane  incremental  diffraction  coefficients  may  be 
expressed  as  follows: 


dE™  (r) 


dz'Ej;, 


^ikr 

4'irr 


sin  9 
sin^0o 


4  sin  —  sm  — 

_ £ _ ^  0 

cos  a  +  cos  <'>0 


(3) 


£JEP0(TM)  sign(-ir  -  <t)o)Eij: 


sin  0 


2  sin  po 


sin^  00  cos  a  -I-  cos  cJiq 


Ho 


dE]^^  (?)  -  dz'Zo 

*  ^ sin  00  4'7Tr 


4  cos  %  cos  ^ 


1 


cos  a  cos  <})o  sin  a 


(4) 


■  [sin  <[)  “  (cos  4)  cos  0  +  cos  a  sin  0  cot  0o)  0] 


(5) 
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and 


‘  iJcf  ^ 

dgPOCTE)  -dz'  sign(-IT  -  <t)o)  Zq  ■■  -| - - - X“ 

•  [sin  <j)4>  ~  (cos  4)  cos  0  —  cos  <J)o  sin  6  cot  Bq)]  (6) 

where  0  <  <[)  ^  217  and 

a  =  cos"'  cos  <})).  (J) 

In  these  expressions,  dE™  (r)  and  dEj°^^^(r)  are  the  incremental  diffraction  coefficients  for  the  total 
and  PO  currents,  respectively,  for  a  half-plane  illuminated  by  a  TM  (E-polarized)  plane  wave,  and  dEp(f) 
and  dE^^^(r)  are  the  corresponding  diffraction  coefficients  for  TE  (H-polarized)  plane-wave  illumina¬ 
tion.  The  half-plane  is  defined  in  terms  of  Cartesian  coordinates  (x,y,z)  by  the  equations  y  =  0,  x  ^  0, 
so  that  the  edge  of  the  half-plane  coincides  with  the  z-axis.  The  origin  of  the  coordinate  system  is  taken 
to  be  at  the  differential  element  (incremental  length)  dz'  of  the  half-plane  edge.  The  directions  of  the 
incident  plane  wave  and  the  point  of  observation  are  given  by  the  spherical  coordinate  angles  (0o,  4>o) 
and  (0,  4>).  respectively  (see  Figure  1).  Ej;;  and  are  the  complex  amplitudes  of  the  electric  and  mag¬ 
netic  components,  respectively,  of  the  illuminating  plane  waves  parallel  to  the  edge  of  the  half-plane  at 
the  origin,  Zq  is  the  free-space  impedance,  and  r  =  |r|,  the  distance  from  the  origin  to  the  field  point. 


OBSERVATION 


Figure  1.  Geometry  of  a  Half-Plane  Illuminated  by  a  Plane  Wave 
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The  half-plane  incremental  diffraction  coefficient  for  the  nonuniform  current  and  TM  plane-wave 
illumination  is  obtained  by  substituting  Eq.  (3)  and  Eq.  (4)  into  Eq.  (2)  and  using  elementary  trigono¬ 
metric  identities  to  simplify  the  resulting  expression  to 


„ikr 


dEn^CTM)  (f)r^dz' 


cos 


2  sin 


<l>o 

2 


I  -ksin 


(8) 


Similarly,  the  half-plane  incremental  diffraction  coefficient  for  the  nonuniform  current  and  TE  plane-wave 
illumination  is  found,  by  substituting  Eq.  (5)  and  Eq.  (6)  into  Eq.  (2)  and  simplifying  trigonometrically, 
to  be 


_  U  Ilf-  t 

Sin  -  ^|cos  -j- 1  -I- sin  -j 


sin  <{)  $  - 


cos  4)  cos  0  -I-  (l  -I-  2  I  cos 


*^0  I  «  I 

Tl  S'" 


sin  0  cot  00 


(9) 


For  arbitrarily  polarized  plane-wave  illumination,  the  half-plane  nonuniform  incremental  diffraction 
coefficient  is  then  obtained  by  summing  Eq.  (8)  and  Eq.  (9): 


dE""  (f() 


>ikrr 


1 


4Trrr  sin  Bqi 


cos 


4>of  I 
2 


+  sinf 


2  sin  0r 
sin  0o< 


Hiz,  / 

signfTT  -  <J)or)  Zq - -  [cos  <})f  cos  0^  -t- 

sin  —  ' 


+  2 1  cos 


<l»0f  I 


2J 


sin  0;  cot  0or ) 


-  sign  (ir  -  <t)of)Zo 


■  a, 

; — jr  sin  <}>r  <}>r 
sin- 


(10) 


with 


a 


cos  * 


/  sin  Q( 
'sin  0or 


cos  <[)< 


(11) 
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We  have  inserted  the  subscripts  "C"  in  Eq.  (10)  and  Eq.  (11),  indicating  that  the  subscripts  are  defined 
with  reference  to  a  local  coordinate  system  with  origin  at  the  incremental  length,  to  clearly  distinguish 
these  locally  defined  quantities  from  those  defined  with  reference  to  the  global  coordinate  system  of  the 
reflector  antenna  introduced  in  the  next  subsection. 


2.2  Transformation  to  the  Global  Coordinate  System  of  the  Reflector 

The  correction  to  the  PO  electric  far  field  of  the  reflector  antenna  is  obtained  as  indicated  above  by 
integrating  Eq.  (10)  around  the  rim  of  the  reflector.  Since  all  quantities  in  Eq.  (10)  are  defined  with  respect 
to  a  local  coordinate  system  with  origin  at  the  differential  element  of  the  reflector  rim,  it  is  first  necessary  to 
transform  Eq.  (10)  to  a  global  coordinate  system  which  we  will  take  here  to  be  that  shown  in  Figure  2. 

The  origin  of  the  global  coordinate  system  is  the  focus  of  the  parabolic  reflector  of  focal  length  F  and 
diameter  D,  and  the  z-axis  is  the  axis  of  the  reflector  directed  towards  the  center  of  the  reflector.  Primed 
coordinates  are  used  to  indicate  integration  points  on  the  reflector  rim,  and  unprimed  coordinates  to 
indicate  the  field  points.  It  is  also  necessary  to  define  the  local  coordinate  system  with  origin  at  the  point 
of  integration  on  the  reflector  rim.  Since  we  are  modelling  the  rim  locally  by  the  half-plane  defined  as 


Figure  2.  Global  Coordinate  System  for  Paraboloidal  Reflector 
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shown  in  Figure  1  (the  incident  plane  wave  in  Figure  1  corresponds  to  the  ray  from  the  feed  at  the  focus 
to  the  rim  point),  we  want  the  local  x-axis  to  be  tangent  to  the  reflector  at  the  rim,  normal  to  the  rim,  and 
directed  inward  from  the  rim;  and  the  local  y-axis  to  be  the  inwardly  directed  normal  to  the  paraboloid  at 
the  rim.  Thus 


icr  =  -sin  yf'  -  cos  y  0' 


« 

and 


ye  = 


-cos  Y  f'  +  sin  y0' 


(12) 


(13) 


with  0',  the  angle  at  the  focus  subtended  by  points  on  the  rim,  given  by 


(14) 


Hence 


te  =  ^  X  9e  =  -4>'  • 


(15) 


We  now  systematically  express  the  locally  defined  quantities  of  Eq.  (10)  and  Eq,  (11)  in  the  global 
coordinate  system.  Starting  with  Qq(,  the  angle  between  the  local  unit  z-vector,  and  the  ray  -r'  from 
the  rim  point  to  the  feed  at  the  focus,  it  is  apparent  that  0of  =  'n-/2  so  that 


cos  0or  =  0 


(16a) 


and 


sin  0of  =  1  . 


(16b) 


Next,  since  <|)or  is  the  angle  between  x;  and  the  projection  of  -  f'  on  the  Xf  yf  -plane  (see  Figure  2),  we  have 


cos  4)or  = 


0f 

-  f '  ■  Xf/sin  0of  =  sin  y 


(17a) 


and 
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0' 

sin  4)of  =  -  f '  •  y</sin  Boe  =  cos  — 


(17b) 


from  which 


1  +  sin  -y 


1/2 


and 


f  1-sin 


01 

2 


2 


l,'2 


(18a) 


(18b) 


Proceeding  to  Q(,  the  angle  between  Zf  and  the  ray  from  the  rim  point  to  the  far-field  point. 


cos  0f  =  -$'  •  ff  . 


But  i(,  the  unit  vector  from  the  rim  point  to  the  far-field  point,  is  parallel  to  R,  the  unit  vector  from  the 
focus  to  the  far  field  point.  Hence 

cos  =  -<!)'  R 

=  (sin  <ti'  X  -  cos  <})'  y)  •  (sin  0  cos  <t>x  +  sin  0  sin  (}>y  cos  0i) 

=  sin  0  sin  (<})'  -  (}>)  (19a) 


and 


sin  0;  =  [1  -  sin^0  sin^(4)'  -  4))]''^ 


(19b) 


where  the  positive  root  must  be  taken  since  0  ^  df  ^  v.  (For  0  =  v/2,  sin  0f  =  |cos(4)'  —  (!))|  .) 
The  cosine  and  sine  of  <})  are  given  by 


cos  <t>(  =  ■  Xf/sin  Q( 


and 


sin  <{);  =  ?;•  ^r/sin  0^  . 
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Again,  replacing  ff  by  R  and  then  switching  to  Cartesian  components  and  simplifying  trigonometrically, 
we  obtain 


cos  <t>f  =  ^-cos  ^  sin  0  cos(<|)'  -  <l))  +  sin  ^  cos  sj/sin  6r 


and 


sin  (})r  = 


A'  6^  / 

-sin  sin  0  cos(4)'  -  4))  -  cos  —  cos  0  /sin 


(20a) 


(20b) 


with  sin  Q(  given  by  Eq.  (19b).  The  angle  a  given  by  Eq.  (11)  is  now  also  expressed  in  the  global 
coordinate  system  since  sin  0or,  sin  6e,  and  cos  <t>r  are  given  by  Eq.  (16b),  Eq.  (19b),  and  Eq.  (20a), 
respectively.  In  particular. 


sin 


a 

2 


( 1  -cos 
\  2 


(21) 


with 


cos  a  =  sin  0r  cos  4)^  . 

The  unit  vector  df  is  transformed  to  the  global  coordinate  system  by  starting  with 

0r  =  cos  0fcos  4>f^  +  cos  0fsin  4>ty<  ■”  sin  dfi(  , 

substituting  global  Cartesian  component  expressions  for  Xf,  ye,  and  ie>  and  simplifying  trigonometrically, 
thereby  obtaining 

0f  =  -  [^cos  0fcos  4>'  (cos  y  COS  4>r  +  sin  y  sin  4>r)  +  sin  0rsin  4>']  * 

-  (cos  0rsin  4>'  (cos  y  cos  4>r  +  sin  y  sin  4>r|  “  sin  O^cos  4>'j  ^ 
sin  4>r  ”  sin  y  cos  4>r 

with  cos  0f,  sin  6e,  cos  4>r,  and  sin  4>r  given  by  Eq.  (19a),  Eq.  (19b),  Eq.  (20a),  and  Eq.  (20b)  respec¬ 
tively.  A  similar  procedure  gives 
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t  /  0^  0^  \ 

=  cos  «t>'  ^cos  -y  sin  «t»f  -  sin  ^  cos  4>f  j  A 

(8'  0*  \ 
cos  ^  sin  -  sin  y  cos  ^tj  y 


(0'  0'  \ 
sin  ^  sin  <1);  +  cos  ^  cos  4><  j  z. 


(23) 


Next,  proceeding  to  e^^Vri  and  replacing  r;  by  r  (see  Figure  2), 


gikrt  eik(R-rR) 

__  =  __  _ 


(24) 


where  R  is  the  distance  from  the  reflector  focus  to  the  far-fleld  oint.  It  is  then  simple  to  find  that 


f'  •  R  =  r'  (sin  0'  sin  0  cos(<|)'  -  4))  +  cos  0'  cos  0] 


(25) 


where 


r'  = 


2F 

1  +  cos  0' 


=  F  sec^ 


and  0'  is  given  by  Eq.  (14).  The  differential  length  dz<  is  given  by 


(26) 


dz(  =  -  r' sin  0'd<j»'  =  -  2F  tan -y  d<t)'  (27) 

and,  referring  to  Eq.  (15),  Ei^<  and  become 

^izf  “  ~  (28a) 


and 


Hjzf  “  Hi,),' 


(28b) 


in  the  global  coordinate  system.  Substituting  explicit  expressions  for  cos  (J>or,  sin  0of,  r(,  dz(,  Ej^^,  and 
Hizf  in  Eq.  (9),  and  noting  from  Eqs.  (17a,  b)  that  0of,  and  hence  also  <i>Qt/2,  lies  in  the  first  quadrant  for 
the  reflector  geometry  we  arc  considering,  gives 
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— .  ikR 

dE"**  (f)  r'  sin  e'd<|)'e“‘'^ 


cos  -:r- 


1 


<l>o<  ■  a 
2 

,'  2  sin  B(  sin  +  ZoHi,),' 

2  2 


E,4.'2 


cin  /K« 


(29) 


where  expressions  in  global  coordinates  have  been  derived  above  for  all  locally  subscripted  quantities. 


2.3  Incremental  Transverse  Far-Field  Components 

Before  integrating  with  respect  to  <{»',  it  is  desirable  to  obtain  explicit  expressions  for  dEg“(f)  and 
dEj“(f),  since,  when  integrated,  these  will  yield  the  0-  and  <|>-components  of  the  nonuniform  current 
electric  far  field  (the  radial  component  of  the  far  field  is,  of  course,  zero).  Thus,  we  let 


CT  =  0#  •  0,  T  =  0r  •  4*.  M-  =  4>t  ■  6i  ■ 


(30) 


Substituting  Eq.  (30)  in  Eq.  (29)  gives 


.ikR 


dEj“  (?)  r'  sin  0'd<l)'e 


'HA'«-ikf’  R 


4>o/  ,  !  a 
cos  -y  +  sm 


.  .  .  .  <bn/  _ _  cos  <t><cos  0ri  „  sin  4>« 

2  sin  Of  sin  —z~  +  ZoHj^,'  - r - 1  <t  —  ZoHi^,’  —  ^  (i. 

2  sin  ^  '  sin  . 


(31) 


and 


aikR 


dEj"  (f)  r'  sin  0'd4.'e-‘«'  R 


4>ot  ,  .  a 
cos  -^  +  sin  — 


4>or  ~  ..  cos  (t){COs  Of 

Eli,'  2  sin  Of  sin  -:r-  +  ZoHi,j,'  - r - 

2  2 


^  „  sin  <1>; 

T  -  ZoHi,(,'  —  ^  V 

sin  2 


To  find  o  and  t,  we  start  with  Eq.  (22),  which  can  be  represented  as 


(32) 


Sf  =  a5i  +  b^  +  ci 


so  that 


cr  =  a  cos  0  cos  ij)  +  b  cos  0  sin  4)  -  c  sin  0 


(33) 
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and 


T  =  —  a  sin  <{)  +  b  cos  4) 


(34) 


Substituting  explicit  expressions  for  a,  b,  and  c  from  Eq.  (22)  into  Eq.  (33)  and  Eq,  (34),  and  making  the 
convenient  change  of  variable 


u  =  <j)'  -  <t) 


(suggested  by  the  <{>'  -  4»  dependence  of  cos  Of,  sin  d(,  cos  (})<,  sin  <();,  a,  and  r'  •  R)  so  that 


cos  —  ^os  <J)  cos  u  —  sin  <|)  sin  u 


and 


sin  4)'  =  sin  4>  cos  u  +  cos  4>  sin  u, 


and  substituting  for  cos  0?,  sin  0f,  cos  4)f,  sin  4)f  from  Eqs.  (19)  and  (20)  leads  to 


cr  =  - 


cos  0  sin  u 


(1  —  sin^  0  sin^  u)*/^ 


and 


T  = 


cos  u 


(1  -sin^  0  sin^  u)'/^ 


Similarly,  by  starting  with  Eq.  (23)  for  4>t,  we  obtain 
cos  u 


H  =  - 


(1  -sin^  0  sin^  u)*/^ 


(35) 


(36a) 


(36b) 


(37) 


(38) 


(39) 
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and 


V  =  — 


cos  9  sin  u 


(1  -sin^  0  sin^ 


(40) 


2.4  Feed  Illumination 

It  is  now  necessary  to  specify  the  form  of  the  feed  illumination  of  the  reflector.  This  is  chosen  to  be 


.tier' 


Ei  =  ^  [a(e')cos  <1>'0'  +  d(0')sin  <{)'<{)'] 


(41a) 


and 


Hi  =  07  [a(0')cos  <1)'4>'  -  d(0')sin  <})’0']  (4lb) 

so  that  the  E-plane  is  the  xz-plane  and  the  H-plane  is  the  yz-plane.  For  an  x-directed  electric  dipole  source 

a(0')  =  cos  0',  d(0')  =  - 1  ,  (42) 

while  for  a  Huygens  source  with  an  x-directed  electric  dipole  and  a  y-directed  magnetic  dipole 

a(0')  =  1  +  cos  0'  ,  d(0')  =  -  (1  -i-  cos  0')  .  (43) 


Thus 

glkr'  pikr' 

Ei,^'  =  -p-  d(0')sin  4)'  =  —p—  d(0')(sin  <})  cos  u  -I-  cos  <J>  sin  u)  (44a) 

and 

gikr'  gikr' 

Hi,),’  =  2^  a(0')cos  <J)'  =  2^  a(0')(cos  <{)  cos  u  -  sin  4)  sin  u)  ,  (44b) 


applying  Eq.  (35).  Substituting  Eq.  (44a)  and  Eq.  (44b)  in  Eq.  (31)  and  Eq.  (32)  and  replacing  d4>'  by  du 
gives 
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ikR  -t 


cos  -^  +  sin  ^ 


2d(0')sin  e^sin  (sin  cos  u  +  cos  «|>  sin  u)<t 


+  a(e') 


cos  6  cos  u  —  sin  6  sin  u 


sin 


(cos  6<  cos  <|>;  (T  -  sin  <j)f  p,) 


(45) 


and 


dE: 


sikR 


cos  +  sin  ^ 


4>o^ 


2d(0')sin  0<sin  (sin  4>  cos  u  +  cos  4)  sin  u)  t 


+  a(0') 


cos  <{>  cos  u  —  sin  4»  sin  u 


sin  I 


(cos  0(  cos  4>f  T  -  sin  4)#  v) 


(46) 


2.5  Integration  of  the  Incremental  Fields 

While  Eq.  (45)  and  Eq.  (46)  can  be  integrated  numerically  as  they  stand,  substantial  simplification 
as  well  as  analytic  clariEcation  can  be  effected  by  eliminating  those  parts  of  the  integrands  that  are  odd 
functions  of  u  and  hence  vanish  when  integrated  with  respect  to  u  from  -it  to  tr.  [Note  from  Eq.  (15) 
that  in  order  to  obtain  the  proper  sign  of  the  far  field  components,  the  integration  must  be  performed 
in  the  negative  u  or  negative  4>'  direction,  since  this  corresponds  to  integration  in  the  positive  Z( 
direction;  see  Eq.  (15).]  It  is  straightforward  to  show  that  all  terms  in  dEj*^  (f)  containing  sin  4>  and  all 
terms  in  dEj"  (f)  containing  cos  4>  are  odd  functions  of  u  so  that 


ES“  (R)  f-«“cos  4>  sin  0'  f  "  du  R) 


1 _ 

4>or  ,  !  ^ 
cos  —  +  sin  - 


2d(0')sin  0r  sin  sin  u  ct  +  a(0')  ^  (cos  0?  cos  4>f  tr  -  sin  4>(  P-) 

2  sin- 


(47) 


and 


EJJ“  (R)  '■■^sin  4>  "^5^  s**'  ®  f 


cos  ^  +  sin  ^ 


4>or 


2d(6')sin  0^  sin  -r—  cos  u  t  —  a  (0')  - —  (cos  0f  cos  4>r  t  —  sin  4>r  v) 

2  sin  2 


(48) 
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Thus  Eg''(R)  and  E^''(R)  vary  as  cos  <j>  and  sin  4>  respectively.  If  the  co-  and  cross-polarized  field  comptonents,  Eco 
and  Ecr.  are  defined  according  to  Ludwig’s  third  definition'^ 


Eco  =  Ee  cos  <})  +  E4,  sin  4)  , 


(49a) 


and 


E„  =  Ee  sin  <|)  -  E^,  cos  4>  , 


(49b) 


then  there  is  no  cross-polarized  field  in  either  the  E-  or  H-planes. 


3.  CALCULATIONS 

In  Section  31,  we  enhance  the  accuracy  of  the  far  field  of  the  paraboloid  reflector  by  numerically 
integrating  the  nonunifbrm  incremental  diflhiction  coefficient  around  the  rim  of  the  reflector  and  adding 
this  correction  or  ‘‘nonunifbrm"  field  to  the  PO  fields  of  the  reflector.  The  E-plane,  H-plane,  and 
4S‘’-piane  amplitude  and  phase  patterns  of  a  20X.  diameter  reflector  with  F/D  =  0.4  are  computed  for 
both  electric-dipole  and  Huygens  sources  feeding  the  reflector.  Plots  of  the  PO  far  fields  and  the  PO  plus 
nonunifbrm  far  fields  are  shown  with  and  without  the  primary  far  fields  of  the  feed  included.  Excellent 
agreement  is  found  with  the  patterns  obtained  firom  the  method  of  moments  (surface  integral  equation) 
solution  of  Bach,  Viskum  and  Frandsen.’-^® 

In  Seaion  3.2,  we  compare  the  amplitude  and  phase  of  the  nonuniform  far  field  calculated  by 
numerical  integration  of  the  nonunifbrm  diffiaction  coefficient  with  that  calculated  by  conventional 
asymptotic  evaluation  of  the  diffraction  integrals.  Expected  good  agreement  is  seen  except  in  the  forward 
and  back  caustic  regions  where  the  asymptotic  expressions  diverge  to  infinity. 

3-1  Far  Fields  of  a  Thin  Metal  Paraboloid  Reflector 

The  field  of  a  paraboloid  reflector  was  calculated  by  integrating  the  nonuniform  current  incremental 
difliraction  coefficient  around  the  reflector  rim  and  adding  the  result  to  the  PO  reflector  field.  The  PO 
field  was  obuined  from  the  following  formulation  due  to  Rusch;'^ 


Ejo  =  ikFcos<},  Vi  ^ 


-cos  9  cos  9’) 


cos  0' 


(a(0')  cos  0  tJo(P)  -  J2(P)]  -  d(0')cos  0  (JoO)  +  J2O)] 


-  2la(0')  sin  0  tan(^)  JiO)j  sin  0'd0'  . 


(50) 
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(51) 


pikR 

—  ikF  sin  <|)  —5- 
K 


gikr'{l  -  cos  0  cos  6') 
l  +  COS  0' 


•  {a(0')  LJoCP)  +  J2(p)l  -  d(0')  Uo(3)  -  J2(3)]}  sin  0'de' 


In  these  equations,  R,  0  and  <J)  are  the  spherical  polar  coordinates  of  the  field  point;  0'  is  now  a  variable 
instead  of  being  fixed  at  the  angle  subtended  by  the  reflector  rim  as  in  Section  2;  0o  is  the  angle  at  the 
focus  between  the  z-axis  and  a  ray  from  the  focus  to  a  point  on  the  rim;  a(9')  and  d(0')  specify  the  feed 
illumination  [cf.  Eqs.  (4la,b)];  Jo  and  J2  are  the  Bessel  functions  of  order  0  and  2,  respectively;  (3  =  kr' 
sin  0  sin  0';  and  r'  =  2F/(1  +  cos  0'). 

To  enable  comparison  with  the  calculations  made  by  Bach  and  Vislcum^’^°  using  a  moment  method, 
the  reflector  diameter  and  focal  length  w'ere  taken  equal  to  20  and  8  wavelengths,  respectively,  and  two 
different  feeds  located  at  the  reflector  focus  were  used;  a  Hertzian  x-directed  electric  dipole,  and  a 
Huygens  source  with  an  x-directed  electric  dipole  and  a  y-directed  magnetic  dipole.  The  co-  and 
cross-polarized  field  components,  E^o  and  E(-r,  were  calculated  from  Eqs.  (49a,  b).  In  the  results 
presented  below,  the  amplitude  and  phase  are  normalized  to  the  amplitude  and  phase,  resp)ectively,  of 
the  sum  of  the  copolarized  PO  and  nonuniform  current  field  in  the  mainbeam  direction  of  0  =  180°. 

First,  we  show  the  results  of  calculations  in  which  the  primarv’  feed  field  is  added  to  the  secondary' 
reflector  field,  as  w'as  done  in  the  calculations  of  Bach  and  Viskum  Figures  3a,  b,  c  show',  respectively, 
the  dipole-fecd  antenna  (co-polar)  amplitude  pattern  for  0  from  180°  to  155°  in  the  E  and  H-planes 
(recall  that  the  cross-polar  field  is  zero  in  these  planes),  and  the  cross-polar  amplitude  pattern  in  the 
4)  =  45°  plane.  Figures  4a,  b,  c  show  the  corresponding  phase  patterns.  Figures  5a,  b,  c  show  the  same 
dipole-feed  antenna  amplitude  patterns  for  0  from  180°  to  0°,  and  Figures  6a,  b,  c  show’  the  correspond¬ 
ing  phase  patterns.  Figures  7a,  b,  c  through  10a,  b,  c  show  the  corresponding  amplitude  and  phase 
patterns  for  the  Huygens-feed  antenna.  In  all  these  patterns,  the  dashed  line  is  the  PO  plus  feed  field, 
and  the  solid  line  is  the  sum  of  the  PO,  nonuniform  current,  and  feed  fields,  so  that  the  difference 
between  the  two  curves  is  a  measure  of  the  importance  of  the  nonuniform  current  field  correction  to  the 
PO  plus  feed  field. 

Comparison  of  our  total  (PO  plus  nonuniform  plus  feed)  field  patterns  with  the  corresponding 
patterns  obtained  by  Bach  and  Viskum  using  a  moment  method,  shown  in  Reference  7  as  well  as  in  a 
more  detailed  and  complete  set  of  curves.^'’  shows  very  close  agreement,  in  general,  throughout.  (Note 
that  in  Bach  and  Viskum  the  main  beam  direction  is  given  by  6  =  0°,  and  an  exp(jwt)  time  dependence 
is  used,  so  that  their  phases  are  the  negative  of  ours.)  As  an  example  we  have  reproduced  in  Figure  1 1 
the  cross-polar  amplitude  patterns  in  the  <t>  =  45°  plane  obtained  by  Bach  and  Viskum^®  for  the 
Huygens-feed  reflector  antenna.  Comparing  Figure  1 1  with  Figure  7c  it  is  seen  that  the  levels  of  the 
sidelobe  peaks  of  our  total  field  pattern  agree  closely  with  those  of  Bach  and  Viskum’s  method  of 
moments  pattern.  Our  null  depths  are  somewhat  deeper  than  theirs,  however.  [This  discrepancy  in 
cross-polar  null  depths  for  the  Huygens  source  may  be  caused  by  the  limited  dynamic  range  of  the 
method  of  moments  program,^”  because  the  cross-polarized  field  in  the  45°  plane  is  the  difference 
between  the  E-  and  H-plane  fields  see  Eq  (498).  and  the  E-  and  H-pIane  fields  for  the  Huygens  source 
are  especially  close  in  the  forward  region  from  0°  to  25°  shown  in  Figures  7c  and  ll.) 
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Figure  3b.  H-Plane  Amplitude  Pattern  of  20\  Dipole-Feed  Paraboloid  Antenna  With  Primary  Field  Included; 
- PO, - PO  +  Nonuniform  Current  Field,  F/D  =  0.4,  0  =  [180°,  155°] 
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Figure  4a.  E-Plane  Phase  Pattern  of  20\  Dipole-Feed  Paraboloid  Antenna  With  Primary  Field  Included; 
- PO, - PO  +  Nonuniform  Current  Field,  F/D  =  0.4,  6  =  [180°,  155°] 
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Figure  4b.  H-Plane  Phase  Pattern  of  20X  Dipole-Feed  Paraboloid  Antenna  With  Primary  Field  Included; 
- PO, - PO  Nonuniform  Current  Field,  F/D  =  0.4,  0  =  [180°,  155°) 
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Figure  4c.  <t)  =  45°  Cross-Polar  Phase  Pattern  of  20X  Dipole-Feed  Paraboloid  Antenna  With  Primary  Field 
Included; - PO, - PO  +  Nonuniform  Current  Field,  F/D  =  0.4,  0  =  [180°,  155°] 
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Figure  5b.  H-Plane  Amplitude  Pattern  of  20\  Dipole-Feed  Paraboloid  Antenna  With  Primary  Field  Included; 
- PO, - PO  +  Nonuniform  Current  Field,  F/D  =  0.4,  0  =  [180°,  0°] 
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Figure  5c.  <j)  =  45°  Cross-Polar  Amplitude  Pattern  of  20X  Dipole-Feed  Paraboloid  Antenna  With  Primary  Field 
Included; - PO, - PO  +  Nonuniform  Current  Field,  F/D  =  0.4,  0  =  [180°,  0°] 
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Figure  6a.  E-Plane  Phase  Pattern  of  20X  Dipole-Feed  Paraboloid  Antenna  With  Primary  Field  Included; 
- PO, - PO  -f-  Nonuniform  Current  Field,  F/D  =  0.4,  0  =  [180°,  0°] 
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Figure  6b.  H-Plane  Phase  Pattern  of  20X  Dipole-Feed  Paraboloid  Antenna  With  Primary  Field  Included; 
- PO, - PO  4-  Nonuniform  Current  Field,  F/D  =  0.4,  0  =  (180°,  0°) 
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Figure  6c.  4>  —  45°  Cross-Polar  Phase  Pattern  of  20X  Dipole-Feed  Paraboloid  Antenna  Wth  Primary  Field 
Included; - PO, - PO  +  Nonuniform  Current  Field,  F/D  =  0.4,  0  =  [180°,  0°] 
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Figure  7a.  E-Planc  Amplitude  Pattern  of  20\  Huygens-Feed  Paraboloid  Antenna  With  Primary  Field  Included; 
- PO, - PO  +  Nonuniform  Current  Field,  F/D  =  0.4,  6  =  [180°,  155°] 
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Figure  7c.  4)  =  45°  Cross-Polar  Amplitude  Pattern  of  20\  Huygens-Feed  Paraboloid  Antenna  With  Primary  Field 
Included; - PO, - PO  +  Nonuniform  Current  Field,  F/D  =  0.4,  6  =  [180°,  155°) 
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Figure  8a.  E-Plane  Phase  Pattern  of  20\  Huygens-Feed  Paraboloid  Antenna  With  Primary  Field  Included; 
- PO, - PO  -t-  Nonuniform  Current  Field,  F/D  =  0.4,  0  =  [180°,  155°] 
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Figure  8c.  <|)  =  45'  Cross-Polar  Phase  Pattern  of  20\  Huygens-Feed  Paraboloid  Antenna  With  Primary  Field 
Included; - PO, - PO  +  Nonuniform  Current  Field,  F/D  =  0.4,  6  =  [180°,  155°] 
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Figure  9a.  E-Plane  Amplitude  Pattern  of  20\  Huygens-Feed  Paraboloid  Antenna  With  Primary  Field  Included; 
- PO, - PO  Nonuniform  Current  Field,  F/D  =  0.4,  0  =  [180°,  0°] 
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Figure  9b.  H-Plane  Amplitude  Pattern  of  20\  Huygens-Feed  Paraboloid  Antenna  With  Primary  Field  Included; 
- PO, - PO  -i-  Nonuniform  Current  Field,  F/D  =  0.4,  0  =  [180°,  0°] 
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Figure  9c.  <\>  =  45°  Cross-Polar  Amplitude  Pattern  of  20X  Huygens-Feed  Paraboloid  Antenna  With  Primary  Field 
Included; - PO, - PO  Nonuniform  Current  Field,  F/D  =  0.4,  0  =  [180°,  0°] 
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Figure  10a.  E-Plane  Phase  Pattern  of  20X  Huygens-Feed  Paraboloid  Antenna  With  Primary  Field  Included; 
- PO, - PO  +  Nonuniform  Current  Field,  F/D  =  0.4,  6  =  [180°,  0°] 
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Figure  10b.  H-Piane  Phase  Pattern  of  20X  Huygens-Feed  Paraboloid  Antenna  With  Primary  Field  Included; 
- PO, - PO  +  Nonuniform  Current  Field,  F/D  =  p.4,  9  =  [180®,  0°] 
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Figure  10c.  <t)  =  45°  Cross-Polar  Phase  Pattern  of  20X  Huygens-Feed  Paraboloid  Antenna  'VTith  Primary  Field 
Included; - PO, - PO  -I-  Nonuniform  Current  Field,  F/D  =  0.4,  0  =  [180°,  0°] 


We  now  make  some  general  comments  on  the  patterns  shown  in  Figures  3  through  10.  Referring  to 
the  dipole-feed  antenna  patterns  of  Figures  3  a,  b,  the  effect  of  the  nonuniform  current  field  in  the  E- 
and  H-planes  is  masked  by  the  primary  dipole-feed  pattern.  In  the  <j)  =  45°  plane  (Figure  3c),  however, 
the  primary  dipole-feed  cross-polar  field  is  close  to  zero  for  0  close  to  180°,  and  the  effect  of  the 
nonuniform  current  field  is  clearly  seen,  resulting  in  a  raising  of  the  nulls.  In  contrast,  for  the  Huygens- 
feed  antenna,  both  the  co-  and  cross-polar  primary  feed  fields  are  close  to  zero  near  the  forward 
direction  and,  referring  to  Figures  7  a,  b,  c,  the  effect  of  the  nonuniform  current  field  may  be  seen  in 
the  E-  and  H-planes  as  well  as  in  the  cross-polar  field  pattern  for  <[)  =  45°.  The  effect  is  most  striking  in 
the  cross-polar  pattern  where  the  nonuniform  current  field  corrects  the  PO  plus  feed  sidelobe  levels 
upwards  by  as  much  as  6  dB. 

Turning  next  to  the  pattern  plot  for  the  full  range  of  0,  it  is  seen  from  Figures  5  a,  b  and  9  a,  b  that 
the  back  lobes  of  the  E-  and  H-plane  patterns  for  both  the  dipole-feed  and  Huygens-feed  antennas  are 
strongly  affected  by  the  nonuniform  current  field.  In  the  E-plane,  the  nonuniform  current  field  raises  the 
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back  lobes,  while  in  the  H-plane  the  back  lobes  are  lowered  from  the  PO  plus  feed  levels.  In  contrast, 
referring  to  Figures  5c  and  9c,  the  back  lobes  of  the  cross-polar  PO  plus  feed  patterns  are  changed  very 
little  by  adding  the  nonuniform  rurrent  field. 

In  the  next  set  of  figures,  we  show  the  results  of  calculations  in  which  the  primary  feed  field  is  not 
added  to  the  secondary  reflector  field.  These  patterns  are  of  more  practical  interest  near  the  forward 
direction,  especially  those  for  the  dipole-feed  antenna,  since  realistic  feeds  are  highly  directive  and  have 
negligibly  small  radiation  in  the  forv'ard  reflector  region.  Figures  12  a,  b,  c  show,  respectively,  the 
dipole-feed  antenna  amplitude  pattern  for  0  fiom  180°  to  155°  in  the  E-  and  H-plane,  and  the  cross-polar 
amplitude  pattern  in  the  <|>  =  45°  plane.  Figures  13  a,  b,  c  show  the  corresponding  phase  patterns. 
Figures  14  a,  b,  c  and  15  a,  b,  c  show  the  same  dipole-feed  antenna  amplitude  and  phase  patterns  for 
the  full  range  of  6.  Figures  16  a,  b,  c  through  19  a,  b,  c  show  the  corresponding  amplitude  and  phase 
patterns  of  the  Huygens-feed  antenna. 
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Figure  11.  <1>  =  45°  Cross-Polar  Amplitude  Pattern  of  20\  Huygens-Feed  Paraboloid  Antenna  With  Primary  Field 

Included; - Spherical  Near-Field  GTD, - Moment  Method, - PO.  F/D  =  0.4,  (After 

Viskum  and  Bach^°) 
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Figure  12b.  H-Planc  Amplitude  Pattern  of  20\  Dipole-Feed  Paraboloid  Antenna,  Primary  Field  Not  Included; 
- PO, - PO  +  Nonuniform  Current  Field,  F/D  =  0.4,  0  =  [180®,  155°] 
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Figure  12c.  (J)  =  45°  Cross-Polar  Amplitude  Pattern  of  20X  Dipole-Feed  Paraboloid  Antenna,  Primary  Field  Not 
Included; - PO, - PO  -f-  Nonuniform  Current  Field,  F/D  =  0.4,  0  =  (180°,  155°] 
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Figure  13a.  E*Piane  Phase  Pattern  of  20\  Dipole-Feed  Paraboloid  Antenna,  Primary  Field  Not  Included; 
- PO, - PO  +  Nonuniform  Current  Field,  F/D  =  0.4,  0  =  [180®,  155°) 
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Figure  13b.  H-Plane  Phase  Pattern  of  20\  Dipole-Feed  Paraboloid  Antenna,  Primary  Field  Not  Included; 
- PO, - PO  -t-  Nonuniform  Current  Field,  F/D  =  0.4,  6  =  (180°,  155°] 
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Figure  13c.  <j)  =  45®  Cross-Polar  Phase  Pattern  of  20X  Dipole-Feed  Paraboloid  Antenna,  Primary  Field  Not 
Included; - PO, - PO  +  Nonuniform  Current  Field,  F/D  =  0.4,  6  =  (180®,  155°] 
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Figure  I4a.  E-Plane  Amplitude  Pattern  of  20\  Dipole-Feed  Paraboloid  Antenna,  Primary  Field  Not  Included; 
- PO, - PO  +  Nonuniform  Current  Field,  F/D  =  0.4,  0  =  [180°,  0°] 
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Figure  14b.  H-Plane  Amplitude  Pattern  of  20\  Dipole-Feed  Paraboloid  Antenna,  Primary  Field  Not  Included 
- PO, - PO  +  Nonuniform  Current  Field,  F/D  =  0.4,  6  =  [180°,  0°] 
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Figure  I4c.  4>  =  45°  Cross-Polar  Amplitude  Pattern  of  20\  Dipole-Feed  Paraboloid  Antenna,  Primary  Field  Not 
Included; - PO, - PO  +  Nonuniform  Current  Field,  F/D  =  0.4,  6  =  [  180°,  0°] 
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Figure  15a.  E-Plane  Phase  Pattern  of  20\  Dipole-Feed  Paraboloid  Antenna,  Primary  Field  Not  Included; 
- PO, - PO  -t-  Nonuniform  Current  Field,  F/D  =  0.4,  6  =  (180",  0°] 
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Figure  15b.  H-Plane  Phase  Pattern  of  20\  Dipole-Feed  Paraboloid  Antenna,  Primary  Field  Not  Included; 
- PO, - PO  -I-  Nonuniform  Current  Field,  F/D  =  0.4,  0  =  [  180®,  0°] 
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Figure  15c.  4)  =  45°  Cross-Polar  Phase  Pattern  of  20\  Dipole-Feed  Paraboloid  Antenna,  Primary  Field  Not 
Included; - PO, - PO  +  Nonuniform  Current  Field,  F/D  =  0.4,  9  =  (180°,  0°] 
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Figure  l6a.  E-Plane  Amplitude  Pattern  of  20X  Huygens-Feed  Paraboloid  Antenna,  Primary  Field  Not  Included; 
- PO, - PO  -I-  Nonuniform  Current  Field,  F/D  =  0.4,  6  =  [180°,  155°] 
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Figure  l6b.  H-Plane  Amplitude  Pattern  of  20\  Huygens-Feed  Paraboloid  Antenna,  Primary  Field  Not  Included; 
- PO, - PO  +  Nonuniform  Current  Field,  F/D  =  0.4,  6  =  [180°,  155°) 
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Figure  l6c.  4)  =  45°  Cross-Polar  Amplitude  Pattern  of  20k  Huygens-Feed  Paraboloid  Antenna,  Primary  Field 
Not  Included; - PO, - PO  +  Nonuniform  Current  Field,  F/D  =  0.4,  0  =  [180°,  155°] 
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Figure  17b.  H-Plane  Phase  Pattern  of  20X  Huygens-Feed  Paraboloid  Antenna,  Primary  Field  Not  Included; 
- PO, - PO  +  Nonuniform  Current  Field,  F/D  =  0.4,  6  =  [180°,  155°] 
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Figure  17c.  4>  =  45°  Cross-Polar  Phase  Pattern  of  20\  Huygens-Feed  Paraboloid  Antenna,  Primary  Field  Not 
Included; - PO, - PO  4-  Nonuniform  Current  Field,  F/D  =  0.4,  0  =  [180°,  155°] 
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Figure  18b.  H-Plane  Amplitude  Pattern  of  20\  Huygens-Feed  Paraboloid  Antenna,  Primary  Field  Not  Included; 
- PO, - PO  +  Nonuniform  Current  Field,  F/D  =  0.4,  0  =  [180“,  0°] 
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Figure  18c.  <|)  =  45°  Cross-Polar  Amplitude  Pattern  of  20X  Huygens-Feed  Paraboloid  Antenna,  Primary  Field 
Not  Included; - PO, - PO  -I-  Nonuniform  Current  Field,  F/D  =  0.4,  0  =  [180°,  0°] 
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Figure  19a.  E-Plane  Phase  Pattern  of  20X  Huygens-Feed  Paraboloid  Antenna,  Primary  Field  Not  Included; 
- PO, - PO  +  Nonuniform  Current  Field,  F/D  =  0.4,  0  =  [  180°,  0°] 
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Figure  19c.  4)  =  45°  Cross-Polar  Phase  Pattern  of  20\  Huygens-Feed  Paraboloid  Antenna,  Primary  Field  Not 
Included; - PO, - PO  -t-  Nonuniform  Current  Field,  F/D  =  0.4,  0  =  [180°,  0°) 


Comparing  Figures  12  a,  b  showing  the  dipole-feed  antenna  E-  and  H-plane  amplitude  patterns 
without  the  primary  feed  field,  with  Figures  3  a,  b  showing  the  patterns  With  the  primary  field,  we  see 
that,  once  the  masking  effect  of  the  primary  field  is  removed,  the  effect  of  the  nonuniform  current  field 
is  principally  to  raise  the  PO  null  depths.  This  is  the  same  effect  that  we  have  already  noted  above  in 
discussing  Figures  7  a,  b  for  the  Huygens-feed  antenna  patterns  where  the  inclusion  of  the  primary 
Huygens-feed  field  does  little  since  it  is  close  to  zero  in  the  forward  reflector  direction.  The  cross-polar 
patterns  without  the  primary  field  dipole  field  shown  in  Figure  12c  are  very  close  to  the  corresponding 
patterns  with  the  feed  field  shown  in  Figure  3c  since,  as  noted  above,  the  primary  cross-polar  dipole 
feed  field  is  close  to  zero  in  directions  close  to  the  mainbeam  of  the  reflector.  The  Huygens-feed 
antenna  patterns  without  the  primary  field  shown  in  Figures  16  a,  b,  c  are  all,  as  expected,  quite  close 
to  the  corresponding  patterns  shown  in  Figures  7  a,  b,  c  with  the  primary  field  included. 


Turning  to  the  dipole-feed  reflector  patterns  without  the  primary  field  for  the  full  range  of  0  shown 
in  Figures  14  a,  b,  c,  we  note  that  in  the  back  directions  (that  is,  0  <  60°),  the  inclusion  of  the 
nonuniform  current  field  has  very  little  effect  on  the  PO  pattern.  This  is  explained  by  the  fact  that  the 
PO  field  in  back  of  the  reflector  is  very  close  to  being  equal  and  opposite  to  the  primary  feed  field  and 
so,  in  the  absence  of  the  primary  field,  masks  the  effect  of  the  nonuniform  current  field.  In  contrast,  the 
effect  of  the  nonuniform  current  field  in  the  back  directions  is  seen  much  more  strongly  in  Figures  5  a, 
b  since  in  these  patterns  the  PO  field  is  largely  cancelled  by  the  primary  field.  A  very  similar  effect  is 
noted  for  the  Huygens-feed  reflector  patterns  without  the  primary  field  shown  in  Figures  18  a,  b  as 
compared  with  the  patterns  with  the  primary  field  shown  in  Figures  9  a,  b.  The  mid-portion  (130°-60°) 
of  the  E-plane  and  the  cross-polar  patterns  of  both  the  dipole-  and  Huygens-feed  antennas  is  strongly 
affected  by  the  nonuniform  current  field. 


3-2  Comparison  With  Nonuniform  Fields  Obtained  by  Asymptotic  Evaluation  of  the  Diffraction 
Integrals 

In  the  previous  section,  we  enhanced  the  accuracy  of  the  PO  far  fields  of  a  20X.  diameter,  8X 
focal-length  reflector  by  adding  the  far  fields  of  the  nonuniform  current  near  the  rim  of  the  reflector. 

These  far  fields  radiated  by  the  nonuniform  current  were  obtained  by  integrating  numerically  the 
incremental  diffraction  coefficients  in  Eq.  (47)  and  Eq.  (48),  which  were  derived  assuming  the  nonuni¬ 
form  currents  near  the  rim  of  the  thin  metal  reflector  are  equal  locally  to  those  of  a  half  plane  tangent 
to  the  reflector  at  the  rim. 

Alternatively,  we  can  evaluate  the  integrals  in  Eq.  (47)  and  Eq.  (48)  asymptotically  by  the  method  of 
stationary  phase  to  obtain  closed-form  approximate  expressions  for  far  fields  radiated  by  the  nonuniform 
current.  Since  the  rim  of  the  reflector  is  a  smooth  circular  edge  many  wavelengths  in  diameter,  the  method 
of  stationary  phase  should  produce  a  very  good  approximation  to  the  integrals  except  near  the  forward  and 
back  directions  of  the  reflector,  that  is,  near  the  caustics  of  the  far-field  rays  emanating  from  the  rim  of  the 
reflector. 

Specifically,  the  two  points  of  stationary  phase  are  simply  the  two  intersections  of  the  rim  of  the 
reflector  with  the  plane  formed  by  the  far-field  direction  and  the  axis  of  the  reflector.  The  contribution 
from  each  of  these  two  stationary  points  of  diffraction  can  be  found  from  the  equations  in  Section  6.2  of 
Kouyoumjian^^  when  the  GTD  diffraction  coefficients  in  Eq.  6.32  of  Kouyoumjian  are  replaced  by  the 
nonuniform  (PTD)  diffraction  coefficients  for  the  half  plane,* 
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Eq.  (52a)  and  Eq.  (52b)  can  also  be  ascertained  by  comparing  D,  and  Dh  in  Eq.  6.32  of  Kouyoumjian 
with  our  Eqs.  (3,8)  and  (5,9),  respectively,  for  0  =  0o  =  tt/Z.  Substitution  of  the  2  -  D  nonuniform 
diffraction  coefficients  Eq.  (52)  into  Eq.  6.30  of  Kouyoumjian  for  the  two  points  of  stationary  phase  on  the 
rim  of  the  reflector  shown  in  Figure  20,  yields  the  following  asymptotic  approximations  for  the  E-  and 
H-plane  fields  radiated  by  the  nonuniform  current: 
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where  all  the  parameters  in  Eq.  (53)  and  Eq.  (54)  have  been  given  previously  except  for  the  fixed  angle 
that  the  tangent  plane  at  the  rim  makes  with  the  z-axis, 
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and  and  E14,  ,  the  incident  magnetic  and  electric  fields  at  <!)'  =  0  and  tt/Z,  respectively,  of  the  electric 
dipole  or  Huygens  feed  illuminating  the  rim  of  the  reflector. 

Figures  21a  and  21b  compare  the  E-plane  amplitude  and  phase  patterns  (radiated  by  the  nonuniform 
current)  computed  from  the  asymptotic  expression,  Eq.  (53),  (dashed  line)  and  from  the  numerical 
integration  of  Eq.  (47)  (solid  line)  for  the  dipole  fed  reflector  20X  in  diameter  with  F/D  equal  to  0.4. 
Figures  22a  and  22b  compare  the  nonuniform  H-plane  patterns  for  the  same  antenna  computed  from 
Eq.  (48)  and  Eq.  (54).  As  expected,  the  agreement  is  excellent  except  in  the  important  caustic  regions 
around  the  forward  and  back  directions  where  the  asymptotic  expressions  diverge  to  infinity.  The  apparent 
discrepancy  in  the  E-plane  amplitude  pattern  (Figure  21a)  near  0  =  55°  is  misleading.  It  is  caused  by  the 
discontinuity  in  the  diffraction  coefficients  along  the  face  of  the  half-plane  when  the  incident  plane  wave  is 
transverse  e|».rtric  [see  Fq.  (SZhl].  This  discontinuity  appears  as  a  large  phase  jump  («  180°)  in  the 
asymptotic  field  Eq.  (53),  yet  as  a  small  amplitude  jump.  The  fields  computed  from  the  numerical  investiga 
tion  in  Eq.  (47),  however,  are  continuous,  so  that  the  180°  phase  jump  near  55°  becomes  a  continuously  yet 
rapidly  varying  field  that  goes  through  a  null  as  shown  by  the  solid  curve  of  Figure  21a.  Fortunately,  the 
efifect  of  the  fields  of  either  the  solid  or  dashed  curves  of  Figures  21a  and  21b  near  55°  on  the  total  field  of 
the  reflector  are  negligible  compared  to  the  PO  fields  near  55°. 
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Figure  21a.  E-Plane  Amplitude  Pattern  of  Nonuniform  Fields  of  20\  Dipole-Feed  Paraboloid  Antenna; 
- Asymptotic,  Eq.  53, - Numerical  Integration,  Eq.  47:  F/D  =  0.4,  6  =  [180°,  0°] 
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Figure  21b.  E-Plane  Phase  Pattern  of  Nonuniform  Fields  of  20\  Dipole-Feed  Paraboloid  Antenna; 
- Asymptotic,  Eq.  53, - Numerical  Integration,  Eq.  47:  F/D  =  0.4,  0  =  [180°,  0°) 
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Figure  22a.  H-Plane  Amplitude  Pattern  of  Nonuniform  Fields  of  20X  Dipole-Feed  Paraboloid  Antenna; 
- Asymptotic,  Eq.  54, - Numerical  Integration,  Eq.  48:  F/D  =  0.4,  0  =  [180°,  0°) 
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Figure  22b.  H-Plane  Phase  Pattern  of  Nonuniform  Fields  of  20X  Dipole-Feed  Paraboloid  Antenna 
- Asymptotic,  Eq.  54, - Numerical  Integration,  Eq.  48:  F/D  =  0.4,  0  =  (180°,  0°) 
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